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Abstract 

The unification of gauge couplings suggests tliat tliere is an underlying (supersymmetric) unifi- 
cation of the strong, electromagnetic and weak interactions. The prediction of the unification scale 
may be the first quantitative indication that this unification may extend to unification with grav- 
ity. We make a precise determination of these predictions for a class of models which extend the 
multiplet structure of the Minimal Supersymmetric Standard Model to include the heavy states 
expected in many Grand Unified and/or superstring theories. We show that there is a strong 
cancellation between the 2-loop and threshold effects. As a result the net effect is smaller than 
previously thought, giving a small increase in both the unification scale and the value of the strong 
coupling at low energies. 
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1 Introduction 



While the Standard Model continues to provide a remarkably consistent description of essentially all 
observed phenomena, many physicists think that there must be a more fundamental theory capable 
of determining many of its features and parameters. To date, the only quantitatively successful 
predictions for physics beyond the Standard Model are the "unification" predictions for its gauge and 
Yukawa couplings and for this reason a great deal of attention has been paid to them. At best, the 
present success of these predictions gives some circumstantial evidence for a stage of supersymmetric 
unification. As the measurement of the couplings of the Standard Model improves, it is possible to 
test such unification predictions with ever higher precision and the hope is this will strengthen the 
case for some underlying unification of the strong, electromagnetic and weak interactions and shed 
light on its specific form. Furthermore, the scale of unification may provide us with the first indication 
of unification with gravity because it is predicted in string theories with a value close to that needed 
for unification. It is the purpose of this paper to test these predictions in detail for a class of models 
which allow for the addition of heavy states as expected in essentially all Grand Unified theories and/or 
superstring theories. 

The "unification" predictions for gauge couplings rest on two premises. The first is that the 
SU{3), SU{2) and U{1) gauge couplings are related at the (large) unification scale and the second 
is that below this scale there is a specific multiplet content. The first premise follows if there is a 
stage of Grand Unification (SU{5), 50(10), Eq etc). However Grand Unification is not essential. In 
a (compactified) string theory the gauge couplings are related even if the gauge group is not unified. 
In this case the relationship between the couplings depends on the Kac Moody level used to construct 
the theory. In the simplest level- 1 theory, the original string constructions yielded the "standard" 
SU(5) relationship between the gauge couplings even though the standard gauge group is just that of 
Standard Mode]|^. Furthermore, in the case of the string theory the scale of gauge coupling unification 
is determined in terms of the only parameter of the theory, the string tension. Fixing this to obtain 
the correct gravitational strength coupling, one finds that the unification scale should be close to the 
Planck scale, i.e of O{10^^GeV). The fact that the unification of the gauge couplings of the Standard 
Model requires a very large scale of unification of 0((1 — 3) x 10^^)GeV may be considered to be the 
first indication that we are on the right track when seeking to extend unification to include gravity 
as in the string theory. However, one should note that the prediction of the scale of unification of the 
Standard Model gauge group couplings applies only to the class of string theory which do not have a 
stage of Grand Unification. 

The second premise needed for the unification predictions is the choice of the low-energy spectrum 
below the unification scale. The usual "minimal" assumption is that this is just the minimal super- 
multiplet spectrum of the minimal supersymmetric Standard Model (MSSM). However this is not the 
most reasonable expectation. 

In any viable Grand Unified theory it is necessary to add additional supermultiplets to obtain the 
required pattern of symmetry breaking. These additional multiplets are vectorlike with respect to 
the Standard Model group, so that they may obtain very large masses of the order of the (Grand) 
unification scale. However, only if they are exactly degenerate with the additional massive gauge 
bosons of the Grand Unified theory, will they leave the "minimal" unification predictions unchanged. 

In string unification, even without Grand Unification, there is a definite multiplet content. We 
noted above that even if the gauge group is just that of that Standard Model, provided it is a level-1 
string theory, the gauge couplings may still be related at the unification scale as in the prototype 
SU{5) GUT. In this case, the breaking of the symmetry is through Wilson lines and one may show 

^Subsequently, alternative constructions were developped which can lead to non-standard U(l) normalisation even in 
level-1 theory - for example see discussion in ref. |l| 
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that the famihes of quarks and leptons belong to complete SU{5) multiplets even though there is no 
residual SU(5) gauge symmetry Q; this is just as is required by observation. In addition, specific 
Wilson line breaking allows for an additional light pair of Higgs doublets as required in the MSSM. 
However, like the conventional GUT, such string theories usually have additional states, vectorlike 
under the Standard Model group, which can acquire masses of the order of the unification scale. In 
cases with a single Abelian Wilson line associated with a freely acting discrete symmetry, these states 
must fill out complete SU{5) representations or come as additional copies of the Higgs sector; no 
other representations are allowed even though there is no Grand Unified group. It is worth remarking 
that the presence of complete degenerate SU{5) multiplets does not affect the unification of gauge 
couplings at one loop and thus, the success of gauge unification applies in this case too. However, 
additional massive Higgs doublets do spoil the unification predictions at one loop and are disfavoured. 

With the minimal spectrum of the MSSM and assuming the SU{5) normalization of the U{1) 
gauge couplings, one may determine the three gauge couplings in terms of the value of the unified 
coupling and the unification scale giving one overall prediction (up to uncertainties due to our lack of 
knowledge of the low energy SUSY spectrum). It is most convenient to express this as a prediction for 
the strong coupling, as this is the least well known of the three couplings. This gives a3{Mz) > 0.126 

including the supersymmetric threshold corrections Q, using as an input the well measured value 
of sin^ 9w and oem- This is in disagreement with the world average value of q;3(M^) = 0.118 it 0.003 
[^, albeit marginally so. In string theories, as we have already mentioned, the gauge unification scale 
is also predicted. In the case of the (weakly coupled) heterotic string the prediction is 

Mstring ~ gstr^n9 X (5.2 X « 3.6 X lO^^GeV (1) 

more than an order of magnitude larger than that predicted by the MSSM. Ways to resolve this 
discrepancy have been recently reviewed by Dienes Of particular interest is a recent suggestion 
of Witten that the (10 dimensional) strongly coupled heterotic string theory (M theory) gives a 
gauge coupling more closely in agreement with the gauge unification value. 

It is clear that the precision measurements of the gauge couplings already allow for a very detailed 
test of the unification predictions and offer the potential of distinguishing between different unified 
theories. Given this, we consider it important to explore in detail the implications for the unification 
predictions of a non-minimal spectrum at the unification scale of the type to be expected in realistic 
extensions of the Standard Model. As we have mentioned, there are likely to be additional vector like 
states with masses of order of the unification scale belonging to complete SU{5) multiplets or copies 
of the MSSM Higgs supermultiplets. Only if these states have exactly the unification scale mass, will 
the usual MSSM predictions apply, but the difference in the predictions if this is not the case applies 
in two loop order and may be expected to be small. We will consider in detail the effects of (5 + 5) 
and (10 + 10) additional SU (5) supermultiplets which have masses less than the unification scale. 
While this does not exhaust the possibilities in Grand Unified theories, as there may be further higher 
dimension representations, the method we develop to obtain precision solutions to the renormalization 
group equations can readily be generalized to such representations. However, in the class of string 
theories discussed above, the only states allowed are indeed just those that fill out the five and ten 
dimensional representations, so the analysis presented here is quite general Furthermore, it is only 
in the context of such string theories that the prediction for the string unification scale may be directly 
compared with the unification scale of the couplings of the Standard Model. 

The organization of the paper is as follows. In the next Section we discuss the mass splitting of the 
additional massive states. This is important because it gives rise to threshold corrections to the gauge 

^ We discard the possibility of additional massive Higgs representations, for they are known to spoil the unification 
predictions at one loop order. 
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couplings. In the case of Grand Unification, the sphtting comes from calculable radiative corrections 
and the resulting threshold corrections are of the same order as the usual two loop corrections. In the 
case of string unification without Grand Unification, the components of an SU{5) multiplet are not 
related by the low energy symmetry and there may be additional threshold effects. 

In Section 3 we develop a method for determining the low energy coupling 03 (M^) and the unifica- 
tion scale, correct to two loop order in the running coupling, ignoring the effect of Yukawa couplings. 
The simplest method is just to integrate numerically the coupled differential equations, technique used 



in and in [10|. However, we found problems in achieving the required accuracy using this technique, 
so we developed a method in which the calculation is done analytically. This has the advantage of 
more clearly exhibiting the important parameters in the calculation. In Section 4 we compare our in- 



termediate results with the form presented by Shifman derived by a completely different method, 
and show that they lead to same answer. In Section 5 we use the analysis to determine analytic for- 
mulae for the unification scale, the strong coupling constant and the intermediate scale masses. Using 
these formulae we determine the numerical predictions and present them in Section 6. The effect of 
additional corrections to the mass of the heavy SU{5) multiplets is considered in Section 7. Finally 
in Section 8 we present a summary of our results and our conclusions. 

2 The Massive spectrum 

We are interested in the class of models with additional vectorlike states I + I where /, / denote 
complete representations of SU{5). The components of complete SU{5) representations transform 
under SU(3) and SU(2) groups as follows: for ip = d'^ : (3,1) , ip = I : (1,2) for the case / is the 5 
dimensional representation of SU{5) and = e'^ : (1,1) ,^p = u'^ : (3,1) ,ip = q : (3,2) for the case / 
is the 10 dimensional representation of SU{5). We note that in string compactifications there are just 
such additional states left massless after compactification . Here we wish to discuss how they obtain 
a (large) mass. Since they come in vectorlike combinations, we can write down SU{3) (8) SU{2) ® U (1) 
invariant mass terms fi^ipip where the mass fi^ can be far above the electroweak breaking scale. This 
term will arise from a stage of spontaneous symmetry breaking through the coupling X^cp^ipip , when 
the SU{3)® SU{2®U{1) invariant scalar field (j)^ acquires a vacuum expectation value (vev). There is 
a very natural explanation for the origin of this vev, because the mass squared of the fields (p^ will be 
driven negative by the radiative corrections involving the coupling A^, in the usual radiative breaking 
mechanism |12| . 

To discuss this possibility in more detail, we consider first the case of a single (j) with Yukawa 
coupling X^cpipip. The renormalization group equations have the form 

^ = (2) 

dirA , ^ / 2 2 2 \ 

= s^y^ [m^ + m^+m^j (4) 

= "2 E Vj,^Mfaj + (ml + m| + (5) 

Here ai = gf/Air , i = {1,2,3}, = A^/ (47r) , b'^ = bi + (ng + 3nio)/2 where bi is the MSSM 
one loop beta function, ?7j_^ are given in the Appendix and Mj are the gaugino masses which evolve 
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as Mi(t) = Mi{0){ai{t)/ai{0)). Also m^, and refer to the soft supersymmetry breaking 
masses of the fields (j), if) and respectively^. The above sum over runs over the component fields 
■0 = {l,dP ,e'^ ,u'^ ^q} for the case one considers both 5 and 10 dimensional complete representations of 
SU{f)) with = i{2n5,3n5,nio,3nio,6nio}. 

Let us discuss how the field cj) actually acquires a vev, and what the expectation for its magnitude is. 
This requires us to solve the renormalization group equations for the soft SUSY breaking masses , 

, m^. It is straightforward to solve these equations numerically in specific cases. To determine the 
general feature, we solve a more simple example, keeping only the QCD gauge coupling and considering 
the case of a single component field ^p which has QCD quantum numbers. Solving eqs leads to 

the result = and 



a{t) 



a(0) 



7m 



m, 



m + 



2r?,-^M|(0) 
^3(2-7^) 



a{ty\ 
a{0)J 



a{t) 
a(0) 



2-7,, 



1 



(6) 



denotes the fixed point value for 
. Note that this equation implies 



where t = ^ln{fi/Mg) and 7^ = s^Y^/b'^, 7^ = 2Y^/b'^ where Yl 
the ratio Y^p/a^ which can be easily computed from eqs.(|2|) and 
that will be driven negative by the radiative corrections, triggering a vacuum expectation value 
(vev) for (j). This is the general behavior for our simplified example with one field ip only. 

In the case that (p has a "flat" potential 0, the scale A at which this happens, corresponds closely 
to the vev that is induced for (p. This then gives a mass to the vectorlike fields I, I coupled to (p. The 
important point to note is that A is determined by the multiplet content and by the allowed couplings 
of the theory; it is not a free parameter. In the case we are considering, it is easy to determine the 
order of magnitude of A using that 7^ oc 1/nio (or l/n^) and 7^ ~ 1. Hence, for large nio (^5) case, 
appropriate to our use of the fixed points in solving these equations, and with M3(0)/m^(0) << 1, 
we get from eq. (^) that: 



lit) 



a{t) 
a(0) 



7m 



m2(0)-2m2(0) 



a{0) 



-7m 



(7) 



The point at which m| is driven negative is clearly dependent on the initial value of mfj,/m^. With 
this of order 1, it will occur at the point where a(t)/a(0) deviates from 1 by 0(1). For the case of 
large number of extra multiplets, this will be very close to the unification scale M^, so A will be close 
to this scale too. 

Although, as we have just shown, it is possible that the intermediate scale breaking is triggered 
close to the unification scale, for a{t)/a{{)) of 0(1), it is also possible for a much smaller intermediate 
scale to be generated, corresponding to a much smaller ratio of a(t)/a(0), its magnitude depending 
sensitively on this ratio. 

Even if these parameters are of the same order, the intermediate scale may be much smaller than 
our initial estimate. If 4> does not have a "flat" potential, then the scale of the vev of (p is not 
determined by the scale at which its mass squared becomes negative. To see this, let us assume the 
potential with leading term c/i'^/M'""^ ^ , where M is the scale of the new physics leading to the effective 
potential, which we take to be the compactification or Planck scale. Then, the vev of (p is proportional 



We use the notation = N5 -\- N^, mo = Nio + N-^ 
*j.e. there are no terms in the potential involving <p. 
^The value of r will be determined by the symmetries of the theory. 
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scale ranging from W^^GeV to W^^GeV for r between 6 and 12. 

To summarize, there is a very natural origin for the generation of the intermediate scale, coming 
from radiative breaking. The magnitude of the intermediate scale is sensitively dependent on the 
details of the theory. 

Let us now consider the expectation for individual intermediate scale masses. 
For the case of Grand Unification, SU{5) or larger, the mass terms must be SU{5) invariant, so 
there is only a single field <j) and is independent of ip. Thus, at the unification scale, all components 
ifj oi a representation / are degenerate. However, below this scale, until the scale M/ =< (p >; 
radiative corrections will split the mass of the components. For the case when one can ignore the 
Yukawa couplings such as in computing the radiative corrections, this splitting is generated by the 
gauge interactions only and is thus easily calculable. The effects of this splitting is included in the 
determination of the gauge unification predictions given in the next Section. The case where Yukawa 
couplings are not negligible is discussed below and in Section 7. 

For the case of string unification, even without Grand Unification, the gauge couplings will be 
unified. Indeed, if the string prediction for the unification scale is relevant, there must not be a stage 
of Grand Unification because the prediction refers to the unification scale for the couplings of the 
gauge group at compactification scale, and this must just be that of the Standard Model. However, 
without Grand Unification, there may be different fields (j)^ responsible for giving mass to the different 
components of the SU{5) representation. Furthermore, there is no reason why the couplings A^ should 
be the same for all ip^s. 

Let us first consider the case that there is only one field (p. In this case the Yukawa coupling X^cpipij) 
leads to a mass \^ < (p > where < > is the vev of 4>. The individual couplings A^ get renormalised 
according to eq(§) causing the spectrum to split. Just as for the SU{b) case, the gauge contributions 
(the first term on the rhs) are included in our analysis below. The case the Yukawa couplings are 
large is discussed in Section 7. This leaves only the corrections which occur if the bare couplings A^ 
are different and the effect of this is also presented in Section 7. 

Finally we consider the possibility when there are different fields <j)^ responsible for the masses of 
the various components tp = {l,d^,e^,u^,q} of the representations of SU{5). Again radiative effects 
will change the couplings A^. However, there is a more important effect, because the radiative breaking 
mechanism giving rise to the <j)^ vev, drives Tn^^ negative at different scales, generating different vevs. 
The terms responsible for this in the renormalization group equations are those proportional to the 
gauge coupling strength of the component field ip. The net effect of this is that the (3, 1) components 
of the 5 dimensional are heavier than the (1,2) components, and the components in the 10 dimensional 
representation are similarly ordered with (3,2) > (3,1) > (1,2) > (1,1). This effect is in the same 
direction as that driven by the gauge radiative corrections to the couplings alone, but generates a 
larger splitting between the components. Again the effect of this may be determined using the results 
of Section 7. 

3 Renormalization group equations 

We consider the two loop renormalization group equations for the gauge couplings, with no top, bottom 
and tau Yukawa interaction. Their general form is given by: 




3 



(8) 



5 



with i = {1,2,3} and where t = ^logQ/Mg] Mg is the unification scale, and bi and bij are, in a 
general notation, the one loop and two loop beta functions respectively, which depend on the multiplet 
content of the theory. Dividing both sides by a? we may rewrite this equation in the form 



-a,; 



-.3 7 

bidt + — y Sdiln aj) + 0{a^ 
4tt . bj 



(9) 



where we have used ajdt 



1 doj 



+ C(a^) to rewrite the second term on the right hand side. We 



integrate this equation from Mg scale down to scale, and replace bi and bij with their appropriate 
values. As discussed in the previous section, the mass splitting of the additional SU{5) multiplets 
added to the MSSM spectrum are generated by radiative corrections. Thus to 0{a^), the threshold 
effects coming from this splitting need only be included at one loop order; the two loop terms may be 
cutoff at the mean mass, /ig of the multiplet. This gives 



a. 



-Si + a~ 



2tt 



with 



1 



1 b - 

4^ bj 



a AM, 



1 

'2^ 



In 



M, 



(10) 



(11) 



In the above equations 6i stands for low energy supersymmetric thresholds; we will express our results 
as a change of the MSSM predictions following from the addition of heavy states and in this the 5i 
cancel; further, ijj stands for the component fields of additional complete SU{5) multiplets we consider 
and runs over {I, d'^, e'^, u'^, q}; their mass is /i^; b'^ and b'^j are the one loop and two loop beta functions 
when none of these states is decoupled]^; their expressions are presented in the Appendix. The above 
formula is correct to two loop order. 

As we have stressed, the one loop threshold effects Q are in fact of the same order as the two loop 
terms. This is made explicit from the one loop renormalization group equations for the running mass 



dt 



(12) 



where r]j^^ are given in the Appendix and we have ignored Yukawa interactions involving the ■0- 
Keeping terms to O(a^) corresponds to keeping only the one loop running approximation for the 
couplings, in the presence of extra-matter: 



dai 
~dt 



Using this we may integrate eq(|l2D to get 



(13) 



In 



^ b', 



a 



(14) 



where we have assumed that all the heavy states get a common mass at unification scale. One sees 
explicitly that one loop threshold effects are, in fact, two loop-like effects. As before we may further 



We assume here that all additional states of a given flavour are degenerate; it is easy to generalise the method to 
the non-degenerate case if desired. 
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simplify this equation to write 



In 



h'- 



OLn 



(15) 



for any ip, since the effect of the radiative mass difference between different flavours generates terms 
of higher order. 

Using eq. (15) to express Q given in eq. (pA|) , we get from eq.(lO) the value of the couplings at Mz 
scale: 



a-\Mz) 



H In 

2n 



47r ^ 6,- 



On 



1 

— y y^,- In 



aj(^<j) 



(16) 



where 



Y 



72 



1^ 

2 6,- 



(17) 



and Ai = 0, A2 = 2, A3 = 3. 

These equations demonstrate an interesting result. At one loop the running of the couplings 
depends only on the number of additional states through the combination n = (ns + 3nio)/2. We see 
that is also true at two loop, even though it is not true for the two loop beta functions, because of the 
threshold effects. Note that this only applies in the case Yukawa couplings are ignored. Note also the 
mild dependence of the functions Yij on n. 

In order to evaluate the effect of the additional heavy states, it is convenient to compare the 
results with heavy states to the minimal MSSM predictions. In the latter, the Renormalization Group 
Equations (RGE) are 



a, 



-\Mz 



-i + -^ln 
27r 



Ml 



^ h or 



(18) 



again ignoring Yukawa couplings. The MSSM variables are labelled with an "o" index to distinguish 
them from the model which considers additional complete SU(5) multiplets. We follow the general 
approach and fix, as an input, the values for a\{M^ and a^iM^) which are known with a good 
accuracy from the experimental values for aemij^z) and sin ^14/ (M^). This fixes the values for a° 
and M° which are then used to compute numerically the value for a\{Mz)- As the MSSM has been 
extensively studied numerically, we can take all these MSSM quantities as known. 

Note that the parameters b\ which appear in eq(|T^) and the parameters bi which appear in eq(^) 
include the constant factors needed for MS — > DK (which are the same in MSSM and our model []), the 
low energy supersymmetry thresholds effects and contributions from nonrenormalizable operators at 
the high scale [O] . The SUSY threshold effects are most conveniently included in terms of the effective 



scale r£,,„ |14] relevant when calculating the strong coupling in terms of the weak and electromagnetic 



" susy 



couplings. Thus, when comparing the prediction with heavy states to the MSSM prediction we find 



ag i(M,) = a°-i(M,) + AaJ ^(M, 



(19) 



'^and will therefore cancel in the calculation below 



7 



where the SUSY threshold effects give a contribution 



susy 



''To 
susy 



and T^"*^ is the new effective scale in MSSM plus the extra heavy states. At one-loop order, the SUSY 
threshold factors are independent of the gauge couplings and thus are the same for the MSSM and 
the MSSM with additional heavy states. At two loop order, the SUSY threshold factors do depend on 
gauge couplings, but the difference between these couplings in the models with heavy states is 0{a) 
and thus the net effect is beyond the two-loop order and can be neglected. The same arguments enable 
us to avoid the explicit form of SUSY thresholds when computing other variables of the model with 
heavy multiplets, like Mg or ^^{^^). 

4 Compatibility with other approaches. 



After completing this analysis our attention was brought to the analysis of Shifman |11| derived from 
a different point of view using holomorphicity arguments and the instanton calculus. Here we wish to 



compare the results of the two approaches, namely the comparison of eq.(16) with Shifman's "master 
formula" We reproduce it below (ignoring the heavy Higgs triplet and 24-plet contributions ^ 

but including the other heavy states' contributions; we also add some generic values, 5j, for the low 
energy Susy thresholds, as we want the values of couplings at scale. We have: 




gen 



In , , £ + 2 In ^ + -ln ^ 



M.ZeR 3 M.ZgL 



+ ^lnJ^ + ^lnJ^ 



3 M.Z^R 3 M,Z, 



dR 



Ma 

+ 2 In ^ 



M,Zh„ 



+ ^ln^ 
2vr fig 



(20) 



+ ttg 



2tt 



In 



+ —\n 



Ma 



2n M,Zh, 



u,d 



Mg 1 



+ ^ln^ 

27T jig 



2 M.ZgL 



1^^ Mg 
2 M.ZiL 



(21) 



03 ^(M^ 



-1 9 1 



M-, 



+ lln. 



2 M.ZdR 



Ma 



(mm.)) 



+ ^ln^ 

2tX fig 



1/3 + 27r ^ 



gen 



M.ZgL 2 M.ZuR 



(22) 



The advantage of this form for the running of gauge couplings is that (above supersymmetric scale, 
when 6i = 0) this form is exact to all orders. However the wave function renormalization coefficients 
Z are only known perturbatively so one is still confined to a given order in perturbation theory 
when testing coupling unification. To two loop order, one must include the values of wave- function 
renormalization coefficients Z in one loop only; also, replacing them with their two loop value would 



These are not included in this work. 
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give a three loop calculation and so on 0. Here we only show that, using their one loop values, we 
get the same two loop form for the running of gauge couplings as we got in equation (|^). The Z 
coefficients are, to one loop, given by: 



Z 



su{j) 



ag 


2C2 

b'. 

J 




2C2 


ag 


2C2 n 

b'. 6,- 
3 ■' 








aj{Mz)_ 




aj{Hq)_ 




aj{Mz) 



(23) 



where j = 1, 2, 3, C2 is the second order Casimir operator for SU{N), and b'j = hj + n and where 
we took into consideration the one loop renormalization of these coefficients. Note that we can neglect 
the splitting of the heavy SU{5) multiplets in the renormalization of the Z factors as they correspond 
to higher than second order corrections in the RGE for the couplingsP], We get that: 



^SU{3) 



^SU{2) 




ag 


2 n 
11 5[ 


ag 


_ai{nq)_ 




[ai{Mz)\ 



2_ 

' 11 



Also we have that 



-'qL 



Z 

ZuR 
ZdR 
ZlL- 



SUiS) 



Z 



SU{2) X ^(/(l) 



1/36 



z 



SUi3) 



z: 



z 



5(7(3) 



X Z, 



4/9 
1/(1) 

1/9 

U{1) 



7 V 7^4 

^SU{2) X ^[/(l) 



Z. 



SU{2) 



X Z 



1/4 

C/(l) 



Z, 



R 



z, 



U(l) 



With these expressions for Z coefficients we get from the "master formula" that: 



(24) 

(25) 
(26) 



-di + Q + — In 
" zvr 



H In 

27r 



fig 



1 o 

I ^ n 



j=i 



25ijXj - 



In 



]2i 



a„ 



aj{M, 



an 



ajinq) 



(27) 



where the bij, bj above are just the values of two loop and one loop beta functions for MSSM. This is 
exactly the formula we previously obtained in eq.(|l6|). 



^Beyond two loop order this form is regularisation scheme dependent JLSf . 

^'^divided by for U(l), hence giving a 3/5 factor, with Y defined as in the Appendix; our definition for Y difi'ers 



byi 

[0 

spe. 

is of higher order. 



from that of and hence the different powers for ^(7(1) in the following equations 

Strictly speaking one should use the bare mass fig in the expressions for Z instead of fj,q, but again, this difference 
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5 Analytical Results 



We now proceed to solve the two loop RGE eqs.(16) analytically to obtain predictions for the model 
with additional SU{5) multiplets. In both the analysis of the minimal MSSM unification and the 
non- minimal case we adopt, we follow the usual analysis method of using the measured weak and 
electromagnetic couplings to determine the strong coupling and the unification scale. Thus, we fix as 
input parameters at the full two loop order ai{Mz) = ai{Mz) and a2{Mz) = a2{Mz) and equal to 
their experimental value In the case of n additional SU{5) heavy multiplets, the mean mass of 
these multiplets introduces a further parameter which we choose to be ag. The output parameters 
of our model will be the unification scale, Mg, the value of alpha strong at the scale, 03 (M^), 
the masses of the additional heavy SU{5) multiplets as well as their common mass fig, at the 
unification scale. Thus, making use of eq.(^), we solve analytically, to two loop order ^ the system 
of equations (|T5|), (|l^) for the unknowns Mg , ij,g , fj,^ , and 03 (Mz) in terms of MSSM variables 
Mg , a^{Mz) and a". This will avoid any numerical evaluation of the system of these equations, as 
well as the inconveniences deriving from possible numerical instabilities. 

To achieve this, we subtract eq.(p!8[) from eq.([l^) for same i and obtain the following relation: 



ar\M,)-ar\M,) 



0-1 / 



a 



-1 



1^9] 


n , 






H In 






27r 





1 

2^ 



Ctr, 



a„ 



oci. 



i3 



\_h_ 

4tt 63 



In 



Q3(M, 



(28) 



Consider now the equations obtained from the above equation for i = 1 and i = 2, and subtract them 
to get: 



h - b2 
2tt 



In 



Mg_ 

MO 



1 b 
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^23 



47r 



will give us the ratio ^ (see later). The eq. (|2^ ) for i 



In 



On 



ajifig) 
a^{M, 



al{M, 



47r ^ 



In 



9. 



(29) 

1 will give us the value of the 



The eq.(|2£ 

variable ^ = l/(27r) ln(Mg//ig) which we introduce instead of fiq as an independent variable and which 
proves useful in solving the system. To see this we make use of the eq.(|l5|) considered for = q: 



- 27r^ - In 



Mr, 



j=i 



OjiHg) 



(30) 



to substitute the ratio jf- in eq. (P8|) . With the above three equations we compute Mg/Mg, fXg/Mg 



(and therefore fig/Mg) and ^ in function of the remaining terms present in the above equations. 
We also observe that ln{a3{Mz)/a2{Mz)) = in two loop order - hence we can drop it - as we 



can see it from the fact that Og ^ 



in one loop order. We also extensively use that 



ln{ag/aj{fig)) = ln(l + ag6j-^) which is correct to two loop order. We therefore end up with expressions 
for Mg/Mg, fig/Mg in terms of ^ and with a nonlinear equation for ^ itself. On this latter equation 

^^The couplings are normalized such as ai — 02 = as at the unification scale. 

^^This means that, in the following , we can discard terms which originate from three loop order in either MSSM or 
in the model with additional multiplets. 
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we again observe that a„ 



o-l 



+ = in one loop order as all log's present are two-loop or 



higher order log's. Hence the ^ dependence of Mg/M°, jig/M" is finally lost, and the three equations 
decouple, as we can always replace ln(l + ag6^-^) by its corresponding value obtained for a ^ computed 
in one loop.PI 

To solve for a^{Mz) we simply multiply eq.(^) for i = 1,2,3 by 62 — &3, &3 — &i and hi — 62 
respectively, add them together and make use of the above observations. In all these processes the 
difference 5i — 5° is dropped as it would bring a higher order contribution, as we mentioned in the 
previous section. 



Finally, the remaining masses of the fields ijj can be computed easily from equation (15). 
After this algebra and consistent neglect of higher order terms, which proves to be essential in 
simplifying the equations, we get the following analytical results: 
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The values of pj , ujj , aj , Tj^^ , with j = {1, 2, 3} and -0 = {h d^, e^, u^, q}, are given by 

(In —39 n 4 n 
" ll2 6['^6[' 216^ 



(33) 
(34) 

(35) 



and finally. 



lln 



-2 n 81 n 2 n 



llvr b[' Utt h'2 7tt b'^ 
7 -3(111 + 13n) 4(22 + n) 



132 6; ' 



28 6' 



21 6' 



6' 



(36) 
(37) 

(38) 



The values of r]j^^ and b'j appearing above are given in the Appendix. 

This gives the full analytical solution, to two loop order, to the system of non-linear equations 
(15),(p!6|),([T7|), which contains three RGE equations and five equations for the running of the masses 
of V^'s. 

The term 1168/231/7r ln(Qg/a°) in the above expressions can be dropped, as it brings in a higher 
order (three loop) correction; however, we keep it, as it gives an idea of the size of three loop corrections 
and makes a better agreement with the numerical approach of this problem. It is also relevant 

We will retain, however, a two loop term in the final results, to get an idea of the importance of three loop 
contributions 
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when considering unification at strong Og. This is equivalent to considering that a~ — 0(g~ + = 
— 1168/231/7r ln(ag/a°) where the right hand side of the latter is a two loop correction 

The solution demonstrates explicitly that the result depends on the number of additional states 
only, via the parameter n = {n^ + 3nio)/2. Also, the values of 03 (M^) as well as the ratio ^ depend, 
to two loop, on the ratio of the couplings at unification scale only. As a check of our calculation, one 
can take the limit a^/a" — > 1 in the above expressions to recover the MSSM case. 

6 Numerical results 

In this section we consider the numerical results given by the expressions of eqs. (|3l|), (|32|), (^3|), (^). 

For n = 1,2 the ratio Mg/M° is rather small, though larger than unity; in the meantime a2,{Mz) 
is always larger than its MSSM value; it cannot became smaller than 03 (M^) because in that case, 
the additional states would be too light 10 — lOOGeF) and these cases are therefore disfavoured. 
Moreover in these cases the two loop approximation breaks down for ag quite close to its MSSM 
value, because the two loop contributions involve log's of very small arguments, which make their 
contribution dominate over the one loop contribution. 

The case re = 3 is special because the step used in deriving eq(|9|) breaks down as it involves 
division by a one loop beta function which is actually 0. This case should be considered separately, 
even though, interesting enough, the limit n = 3 in our final results exists and is finite, predicting a 
larger unification scale, up to a factor of ~ 10 for unification at strong coupling together with a larger 
0:3 (M^). However the higher order corrections are expected to be very large in this case too so these 
results should be viewed with caution. This case has been analysed in detail ||l^, but ignoring the 
gauge corrections to thresholds which, as we have seen, affects the result considerably. 

For n = 4 or larger the perturbation expansion is well behaved for intermediate values of the 
coupling. The results for the unification scale and the strong coupling are shown in Figs 1 and 2. One 
may see that the unification scale is increased, but by a factor < 5 for Ug < 0.45. The values of 03(^/2) 
always increases and are larger than in MSSM. This is in disagreement with the analysis of |jl^ due, 
we think, to errors in the numerical intergration used. In fact it can be shown analytically that the 
partial derivative of 03 with respect to the ratio ctg/ctg is always positive for any fixed re > 4 with 
otg > otg (see its expression given in the Appendix). Therefore, in these cases, the strong coupling is 
larger then in MSSM. 

The dependence of the masses of the additional states on the unified coupling is shown in Figs. 3, 
4 and 5. One may see the masses increase as we increase the number of states re, while keeping Og 
fixed. Note also that for a given n, there are two values of Og which predict the same value of /ig, 
which shows the advantage of setting Ug as our input parameter instead of fig , avoiding the ambiguity 
of the methods having ^g as an input parameter and Og as an output. This is particularly relevant 
when we consider the predictions for large values of the unified coupling. 

To summarize, in contrast to previous analyses, we have found that the addition of heavy states 
systematically increase 03 (M^) taking it further from the current experimental value. The effect of 
the additional states does increase the unification scale, but by a factor less than ~ 5 even for the 
largest couplings which may still be in the perturbative domain. Since for unified couplings of 0(0.3) 
the string prediction of eq(|^) increases by a factor 2.7 there is only a small improvement with the 
weakly coupled heterotic string prediction. Note that the analysis has neglected the Yukawa couplings 
of the third generation which may be important; their effect is investigated elsewhere [ p!7[ . 

^^See Appendix for a full two loop equation 
i^We use a'^{M^) = 0.125 
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7 Massive threshold corrections 



The analytic results we obtained in Sections 5, 6 refer to the case when the heavy multiplets have a 
common bare mass, fig at the unification scale and are split only by gauge interactions. We found that 
the threshold corrections introduced by this splitting largely cancel the two-loop corrections in the 
beta functions leading to the form of eq(|l^) or equivalently eq(^). However Yukawa couplings will 
also cause the massive spectrum to split introducing further threshold corrections. Fortunately the 



form derived by Shifman, eq(22) shows that the cancellation persists for Yukawa corrections as well 
and that the heavy states contribute with their bare mass only; radiative corrections to the masses 
should not be included if the appropriate two-loop beta function is used. Thus the only effects of 
Yukawa couplings are contained in the Z factors for the light fields. This is important because, unlike 
the threshold effects, these corrections are not proportional to n and so we do not expect such effects 
to be large. There remains the corrections which may occur in non-Grand-Unified models due to the 
possible splitting of the bare Yukawa couplings introduced in Section 2. In this case the correction 
factors to our previous results are given by 
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(39) 



where the label "old" used in this section refers to our previously obtained results of Sections 5 and 
6 and //^ stands for the bare mass for ip field. Also: 



a^\M,)-a^\M,)M = ^ln 



, 3nio . 

H m 

287r 



o5,,o2 



(40) 



for the change in the strong coupling. The general effect of this uncertainty requires a specific model 
for the bare couplings. However a general feature emerges. For the case of additional 5s an increase 
in the unification scale is always accompanied by an increase in the strong coupling. For the case of 
additional 10s it is possible to avoid this, but only through a very special choice of the bare masses 
which, lacking a specefic model, looks quite unnatural. Further one may see that the change in the 
unification scale is quite small for reasonable numbers of additional multiplets unless the bare masses 
are split by many orders of magnitude. In string theories the differences between the bare Yukawa 
couplings are typically of 0(1) for the allowed couplings so it seems unlikely such bare threshold 
corrections will give a large change to the unificaion scale. 



8 Summary 

Motivated by the structure of many Grand Unified and Superstring theories, we have considered 
the implications of extending the structure of the MSSM by the addition of extra copies of massive 
vectorlike (/ -|- /) representations. These we take to fill out complete SU{5) multiplets, even though 
there may be no stage of Grand Unification. While such representations do not affect the one loop 
unification predictions, they can have a substantial effect through their threshold effects and at the two- 
loop level. We computed their effects on the unification predictions for the running gauge couplings for 
the case the additional representations have common bare mass. We found that they always increase 
the value of as (ignoring Yukawa couplings effects) and thus, cannot explain the (marginal) discrepancy 
between the predicted value as{Mz) > 0.125 and the observed value a^[Mz) = 0.118 ± 0.003. 

We also found that the unification scale is systematically increased by the addition of such heavy 
states. The increase depends on the value of the coupling at unification and is largest for the case 
this is approaching the strong coupling limit. For this case the unification scale increases by a factor 



13 



varying (function of n) between 1 and 5 for ag=Q.3, at the limit of the perturbative analysis. The 
larger enhancement, allowing for the change in the string prediction at larger coupling, is only slightly 
closer to the value predicted in the (weakly coupled) heterotic string theory. 
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10 Appendix 

The one loop beta function in the presence of complete SU (5) multiplets is given by: 



b' = 



( f + n 



V 




(41) 



where n represents the linear combination n = (n5 + 3nio)/2, with ns = N^+N-^ and nio = -^io+-^io- 
The one loop MSSM beta function is given by 6j = 6'j(n = 0). 

The two loop beta function in the presence of the complete SU{5) multiplets is given by: 



/ 199 , J7_r) -I- 



23„ 27 , 9 „ I 3 „ 
nio -5- + io»^5 + Tn'^10 



30 '"f ' 10'"^^ 5 ' 10 '"f ' 10 

I + ^"-5 + i7t"-io 25 + |n5 + ^nio 
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^ + + fnio 



9 + 3n 
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24 + 8nio 
14+^n5 + 17nio / 



(42) 



The MSSM two loop beta function is given by bij = 6^j(n5 = 0, nio = 0)- 

The values of the coefficients Pj^^p used in text, are given by the following matrix with the symmetry 
group index j as a line index j = {1, 2, 3}, and with the type of the field ip as a column index running 
over the set {I, (F^ e'^, u"^, q}, in this order. 



In the same notation we have for the matrix r/j_^ 




(43) 



( ^ 
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V 



2^ 

15 





6 

5 






_8_ 

15 





(44) 



where Cji^p) is the quadratic Casimir operator, given by 3y^/5, 3/4, 4/3 for C/(l), SU{2), SU{3) 
respectively. 

With these definitions, we get. 



V'={«,d=,e=,u=,q} 



(45) 



14 



with Ai = , A2 = 2 , A3 = 3. 

Throughout the text we used SU{5) normalization for gauge couplings, t = l/(27r) ln(sca/e/Mg), 
a, = gf/{A7r),Y^ = Xl/{47r). 

In two loop, the equation for used in Section 4 is 
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'23l7r 
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with A J given by 



A. 



13n 291n -24n 



13206; tt' 566(,7r' Ib'^ir 



(46) 



(47) 



The derivative of 03(^^2) (eq.(|3^)), with respect to 2; = ag/a" is positive, given that z > 1 and n > 4 
as it can be seen from its expression given below: 

^asiM^) = 3 [ 320 + 240 (-4 + n) + 60 (-4 + nf + 5 (-4 + nf + 2720 {-I + z) 

+1600 (-4 + n) (-1 + z) + 290 (-4 + nf (-1 + z) + 15 (-4 + nf (-1 + z) 
+5716 (-1 + zf + 2789 (-4 + n) (-1 + zf + 400 (-4 + nf (-1 + zf 
+15 (-4 + nf (-1 + zf + 496 (-1 + zf + 1429 (-4 + n) (-1 + zf 
+170 (-4 + nf (-1 + z)^ + 5 (-4 + nf (-1 + zf] x 
{147r [n + 5 (-1 + z) + (-4 + n) (-1 + z)]x 
[20 + 5 (-4 + n) + 53 (~1 + z) + 5 (-4 + n) (-1 + z)] x 

z[-l + n + (-4 + n) (-1 + z) + z]}'^ al{M^) (48) 
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Figure captions 

Figure 1. The ratio Mg/M° plotted in function of the ratio ag/a" for different values of n. 

Figure 2. The values of a^^M^) plotted in function of the ratio ag/a° for different values of n. 

Figure 3. The values of logiQ{fiec / M") plotted in function of the ratio ag/a" for different values of n. 

Figure 4. The values of logio(Mg/-^g) plotted in function of the ratio ag/ag for different values of n. 

Figure 5. The values of the splitting logio(/^g//^e=) plotted in function of the ratio ag/a" for different 
values of n. 
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